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1 Introduction 

This paper aims at presenting a few models of quantum dynamics whose description involves 
the analysis of random unitary matrices for which dynamical localization has been proven to 
hold. Some models come from physical approximations leading to effective descriptions of the 
dynamics of certain random systems that are popular in condensed matter theoretical physics, 
whereas others find their roots in more abstract considerations and generalizations. Although 
they may differ in detail, the operators describing the models all have in common the following key 
features on which their analysis relies heavily: their dynamics is generated by unitary operators 
on an infinite dimensional underlying Hilbert space which have a band structure when expressed 
as matrices in a certain basis and the randomness of the models lies in phases of the matrix 
elements. 

The focus of this note is put on the description of the models and of the localization results 
available for them. The methods and tools at work in the detailed proofs of these results are only 
briefly presented, with an emphasis on the similarity with the methods used in the self-adjoint 
case. A detailed account of such proofs can be found in the paper [18J to which the reader is 
referred for more about technical issues. 

The paper starts with a model of electronic dynamics that we call the magnetic ring model and 
which, in a certain sense, is the root of the other models that follow. The next section makes the 
connection between the evolution operator of the magnetic ring model and the CMV matrices, 
which play a major role in the theory of orthogonal polynomials with respect to a measure 
on the unit circle. Then we introduce the unitary Anderson models as natural (i- dimensional 
generalizations based on the structure of the evolution operator stemming from the magnetic 
ring model, and on its similarity with the well known discrete Anderson model. A final section is 
devoted to a model of one-dimensional quantum walk in a random environment, another rather 
popular topic of study in theoretical physics and computer science. For all these models, we 
state dynamical localization results which are based on the methods that we describe in the last 
section section of this paper. 
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2 Magnetic Ring Model 



Consider an electron in a metallic ring threaded by a time dependent magnetic flux at the center 
of the ring. Further assume the flux grows linearly with time. According to Maxwell's laws, 
the flux induces a constant electric force tangent to the ring. Hence the electron is submitted 
to the field force induced by the periodic metallic background plus the constant force induced 
by the magnetic flux. A natural question addressed in [271 IHl E] consists in asking whether, 
asymptotically in time, the electron will acquire an unbounded energy due to the constant force 
it feels or if the defects of the metallic structure of the ring can prevent the growth in energy. 

In order to tackle the problem, the following approximations and regime are considered: the 
curvature and width of the ring are neglected and the strength of the constant force is small. 
This leads to an effective one dimensional periodic model in the angular variable, x G [0, 2tt), see 
figure [1} The corresponding Hamiltonian takes the form 

Magn. flux A $(t) ~ at 




Period, pot. Tangent force a 

Figure 1: The magnetic ring model 



H(t) = {-td,-atf + Vpix), on L\{0,27i]), (1) 

with periodic boundary conditions, where the parameter a is assumed to be small and Vp is 
real valued. Note that the variable at plays the role of the quasi-momentum for the periodic 
Schrodinger operator with potential Vp extended to M by periodicity. Therefore the spectrum of 
H(t) is given by the corresponding band functions and is periodic in t. Moreover, the 

effective Hamiltonian being slowly varying in time for a ^ 1, the adiabatic theorem of quantum 
mechanics states that an initial condition proportional to an eigenstate of H{0) will give rise 
at any later time to a solution which belongs to the corresponding eigenspace of H{t) obtained 
by continuity in time, to leading order in a, provided the eigenvalues Ek{t) are simple for all t. 
Therefore, over a period, such an initial state only changes by a phase which depends on the 
potential Vp. In order to describe energy growth, it is necessary to allow transitions between the 
(instantaneous) energy levels of the model. For a quantitative approach, one makes use of the 
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Landau- Zener formula which says that the amphtude of non adiabatic transitions between two 
levels is appreciable only when the gap between the levels is small, actually of order -y/a, so that 
the levels experience an avoided crossing. Now, considering that typically over one period in t 
each level becomes close to the level immediately above and immediately below only once and 
at different times, (except for the ground state), see figure [21 the effective evolution operator is 




Figure 2: The energy levels of H{t) and the different transitions considered 



constructed on the basis of the considerations above as follows: over the first half period, the 
two levels with indices 2k and 2k + 1, k > 0, exhibit one avoided crossing during that time span 
and evolve independently of the others, according to some scattering process. Over the next half 
period, the same scenario takes place, except that the set of independent levels involved in an 
avoided crossing carry indices 2k — 1, 2k, (except for the ground state). For a given set of two 
levels exhibiting an avoided crossing, with indices k — 1, k, with k > 1, the scattering process is 
encoded in a general 2x2 unitary matrix 



Sk 



itkC 



«7fc 



with ak,jk,Ok e [0,27r), 
and Tfc, tk e [0, 1], s.t. + tl = 1. 



(2) 



The coefficient tk gives the Landau-Zener transition amplitude associated with the avoided cross- 
ing and depends only on the minimum gap displayed by the band functions and their local be- 
havior there. The phases depend in a more complicated way on the global behaviour of the band 
functions. When k = 0, Sq is replaced by a phase, sq. In principle, once Vp is given, all param- 
eters of Sk can be computed, within the framework and approximations adopted, see [2Zl El E]- 
Altogether, the effective evolution operator over one period, also called monodromy operator, 
takes the following matrix form on /^(N) in an orthonormal basis of eigenstates of H{0) 



(Si 



U = UoUe, where Uo 



( So 



So 



s. 



\ 



■ J 



(3) 
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We shall denote by {efcjfcgN the chosen basis such that H{0)ek = Ek{0)ek, k eN. Note that the 
2x2 blocks in f/g are shifted by one with respect to those of Uq along the diagonal, and that Sq 
represents a 1 x 1 block. Without expliciting the elements, we have the structure 



U 



\ 



* 
* 



V 



(4) 



Let us note here that not all phases appearing in the matrix U play a significant role. Indeed, it 
is shown in [8] that a suitable change of phases of the basis vectors amounts to setting all phases 
{7fc}fcGN to zero. 

On the basis of the arguments leading to the operator U describing the evolution over one 
period, the large time behavior of the electrons in the ring threaded by a linear magnetic flux is 
encoded in the properties of the discrete dynamics generated by U. This is the starting point of 
the analysis and we shall not attempt to justify rigorously any of the arguments outlined above. 

Coming back to the original motivation, we assume that the periodic potential Vp contains a 
random component due to the impurities in the metallic ring. Then all matrices Sk are random 
and, in turn, the monodromy operator becomes a random unitary operator with band structure 
that we denote by U^^. The subscript u indicates some configuration of the random parameters. 
We will specify below the way the monodromy operator depends on the randomness. Since we 
are working in an energy eigenbasis {ek}kGN, the question asked at the beginning of this section 
can be cast into the following form. 



Question: Let ip G /^(N) be normalized with compact support. 



z.e. 



(cklv) =0,ifk>R, 



for some i? > 0, so that its energy is bounded above by -Eij(O). For a typical configuration of 
impurities u, does the random vector at time n , U^(f, travel to high energy states or spread 
significantly over high energy states of the basis e^. A; > as n — t- oo ? Or does the vector U^ip 
remain close to a finite dimensional subspace spanned by basis vectors Ck with k < p, uniformly 
in n ? A related but not equivalent question is: for a typical configuration u, does the spectrum 
of the operator U;^ contain a continuous component or is it pure point? 

We will be able to provide a quantitative answer this question, for certain choices of deter- 
ministic and random parameters in the model. We shall refrain from stating results in their full 
generality, referring the interested reader to the references provided for more details. Several 
such choices are studied in [5l El [HI |T7l [291 HI]--- We will only discuss one of them which, on 
the one hand, is rich enough for our purpose, and, on the other hand, was actually proposed to 
study the physical model [SI [3] • This model is defined as follows: 
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We assume the transition amplitudes between neighboring levels are deterministic and all 
take the same value, whereas the phases of the scattering matrices are random. This hypothesis 
is certainly a simplification but it also makes the problem more interesting, in the sense that 
transitions to higher and lower energy levels are equally probable, independently of the energy. 
Therefore the random phases through their interferences play the key role. See [5l [8] for discus- 
sion of cases with variable transition amplitudes. 



Assumption A: 

The coefficents {tk, r^) in ^ all take the same value (t, r) G (0, 1)^, for all k > 0. 

We also exclude the trivial case t = such that U^^ is diagonal, and r = such that the 



absolutely continuous spectrum of U^) coincides with the unit circle S, see Remark 2^ below and 

Next, we assume the randomness enters the operator [/^ through phases which are i.i.d. on 
the unit circle. We formalize this as follows. Let [Q, J-", P) be a probability space, where Q is 
identified with {T^}, T = ]R/27rZ being the torus, and P = (g)fceNPA;, where Fk = for any k gN 
and u is a. fixed probability measure on T, and J-" the a-algebra generated by the cylinders. We 
define a set of random variables on {Q, J-", P) by 

Ok-.n^T, s.t. O'^ = tuk, keN. (5) 

These random variables {6'fc}fcgN are thus i.i.d. on T. 
Assumption B: 

Let Di_j = diag {e"*^^} in the basis {ek}k£N, where the 6'j^'s are given in ([5|. Suppose (iz/(r) = 
T{9)d9, where < r G L°°([0, 27r)). 

Under Assumptions A and B, we consider operators f4j of the form 



U^ = D^S, withD^= diag{e-*^^} 



(6) 



and 



/ r 



S 



rt 

rt 
-t^ 



-t^ 
—rt 

—tr 



\ 



rt 
rt 



-t' 
—rt 

J2, 



(7) 



\ -t^ -tr ■•./ 

In the case where all the (relevant) phases in the scattering matrices are i.i.d. and uniform 
on the unit circle, it can be shown that t^j takes the form ([6]) with a uniform density r, see 0. 
This special case is argued to be physically relevant in [9], but the result below holds for any 
density r satisfying assumption B. Note that the operator S is obtained by formula ([3]) with 
blocks Sk of the form 



2fc+l 



r t 
—t r 



2(fc+l) 



-t 

r 



\fk G N, and Sq = 1. 
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Theorem 2.1 JJFJ Consider defined in under assumptions A and B. Let t G (0,1) he 
arbitrary and denote by E the expectation over u. Then there exist a > 0, C < oo such that 



E 



sup|(ej|t/>fc)| 



Consequently, for any p > 0, we have 



sup ||X^?7"v9|| < oo almost surely, 



(9) 



(10) 



where the operator X is defined by Xck = kck, for all A; G N. Moreover, the spectrum of U^j is 
pure point: 

<^{U^) = o'pp{Ui^) almost surely (11) 
with exponentially decaying eigenf unctions. 

The previous statement is a dynamical localization result in energy space. Further assuming 
that Ek{0) < C'kP, as A; — ?• oo for some C ,p < oo, it shows that the energy of the electron in the 
disordered metallic ring does not grow unboundedly with time, despite the constant force acting 
on it. Also, the probability to find the electron in high energy states, i.e. with high quantum 
number number, decays faster than any inverse power of the quantum number. Note however, 
that there are different circumstances where the spectrum of U may be pure point but the energy 
can grow in time, [H]. 

Remark 2.2 It is often technically simpler to consider that the operator C4; acts on rather 
than on This means that all indices k are considered as elements of Z instead of N, that 

fl = {T^}, P = <^kez^, and that we deal with unitary operators of the form 



U... 



and 



S 



D^S, with = diag {e *^^} 



(12) 



v 



rt 








^2 


—rt 






rt 




rt 


-t^ 


-t^ 


—tr 




—rt 






rt 


1^2 






-t" 


—tr 



(13) 



where the translation along the diagonal is fixed by {e2k-2\Se2k) = —t^, A; G Z. 

In particular, on /^(Z), one sees rightaway that if r = 0, f/^ is unitarily equivalent to a direct 
sum of two shifts. Hence it has purely absolutely continuous spectrum given by §. Since one can 
pass from Uui defined on /^(Z) to two copies of the monodromy operator defined on /^(N) by a 
finite rank operator, this shows that aa.dUi^) = § in either case. 



Theorem 2.1 applies to this setting as well, mutatis mutandis., as discussed in [HI [T7] . 
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3 Orthogonal Polynomials on the Unit Circle 



Before we turn to other generalizations of this model, we briefly mention in this section that 
unitary operators with a band structure of the form ^ appear naturally in the theory of or- 
thogonal polynomials on the unit circle. For a detailed account of this topic, we refer to the 
monograph [34j. Given an infinitely supported probability measure dfi on S, such polynomials 
are determined via the recursion 



with $t 



z'^^kil/z), $0 = 1, 



(14) 



by a sequence of complex valued coefficients {ak}k<^n, such that \ak\ < 1, called Verblunsky 
coefficients, which also characterize the measure dfi, see [31]. This latter relation is encoded in 
a five diagonal unitary matrix C on /^(N) representing multiplication by z G §: the measure 
dfi arises as the spectral measure ;u(A) = (eo|-E(A)eo) of the cychc vector cq of C, where dE 
denotes the spectral family of C. This matrix is the equivalent of the Jacobi matrix in the case 
of orthogonal polynomials with respect to a measure on the real axis, and it is called the CMV 
matrix, after [T2]. 

Writing the Verblunsky coefficients as 



ak = re^^'' , and setting 
the corresponding CMV matrix reads 



1 - 

' A:' 



/c = 0,1, 



(15) 



C 



to 



\ 



_7.g7.^g«{»?o-m) 
tlt2 



toil 
-rir2e 



t2h 



(16) 



which is a special case of ([s]), see e.g. [T^. In the same way as tri-diagonal Jacobi matrices can 
be seen as paradigms for self-adjoint operators, the result of [12] shows that five-diagonal unitary 
matrices (|3]) are paradigms of unitary operators. This gives a model independent motivation for 
the study of such operators. 

Comparing with (|6]), it was noted in [17] that if the Verblunsky coefficients all have the same 
modulus and if their phases rj^ = 0^ + 9k-i + ■ ■ ■ + 6*0, then C is unitarily equivalent to —U. 
Therefore, assuming the 9^ are i.i.d.. Theorem 1 then directly yields the 



Corollary 3.1 JT^ E 

Let aA:(i^)fegNg be random Verblunsky coefficients of the form 



k 



0,1,2, 



ak[u)j = re-"""' , < r < 1, 
whose phases are distributed on T according to 

Vki^) ^ du * du * ■ ■ ■ * du , (k + 1 convolutions) 



(17) 



:i8) 
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where dv satisfies assumption B. Then, the random measure djj,^ on § with respect to which the 
corresponding random polynomials ^k,ui ore orthogonal is almost surely pure point. Moreover, 
both and (7^ hold. 



Remark 3.2 Other dynamical localization results for random polynomials on the unit circle 
are proven for independent Verblunsky coefficients, 13^ [3B[ The results of [3^ and 

require rotational invariance of the distribution of the Verblunsky coefficients in the unit disk. By 
contrast, the corollary above holds for strongly correlated random Verblunsky coefficients. 



4 Unitary Anderson Models 



When the unitary operator f/^^ = Di^S is considered on /^(Z) according to Remark 2.2, the 



similarity with the well known (self-adjoint) one-dimensional discrete Anderson model is evident: 



The 2-translation invariant unitary operator S given in (13) plays the role of the translation 
invariant discrete Laplacian A and the diagonal random matrix Di^ is similar to the diagonal 
random potential operator V^. The sum —A + K; is replaced by the product D^S, since we 
deal with unitary operators. Although [/^ 7^ g-«(A+VL)^ ^j-^jg operator can be viewed as an 
effective generator of a discrete dynamics of a particle on the one dimensional lattice. In that 



case, Theorem |2.1| can be interpreted as dynamical localization result in a one dimensional 
configuration lattice, which begs to be generalized to arbitrary dimension. Such a generalization 
was proposed in [20] which we now describe. 



To define the multidimensional version of the unitary equivalent of the Laplacian, we view 
f'il/) as ®^^i/^(Z) and define the canonical basis vectors e^, for k & 1/ hj Ck — ® ... ® e^^. 



Making explicit the dependence in t in 5 = S{t) from (13), we define Sd{t) by 



Sd{t) = ®,ti5(t). (19) 

We denote by | ■ | the maximum norm on Z'^. Using this norm it is easy to see that Sd{t) inherits 
the band structure of S(t) so that 

{ek\Sdit)ei) = if|A;-/|>2. (20) 

Due to the tensor product structure, the spectrum of Sdit) is obtained from that of Sit), which 
can be determined by using Fourier transform. We get 

a{Sdit)) = {e'^ : ^ E [-d\o, dXo]}, where Aq = arccos(l - 2t^). (21) 

The random operator D^^ keeps the same form in the canonical basis, = diag {e~*^fc}, 
with the understanding that {0'^}kez'' are i.i.d. on T, with distribution du. 

The operator 

= D^Sdit) defined on fiZ'^) (22) 
is called the generator of the unitary Anderson model. 
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In that framework, Theorem 2A_ is a unitary version of the statement that dynamical lo- 
cahzation holds true for any disorder strength in one dimension for the Anderson model with 
absolutely continuous distribution of potential. As is well known, localization results for the 
Anderson model in two and higher dimensions are only available in certain asymptotic regimes 
of the parameters, typically large disorder, or in certain subsets of the spectrum, the band edges. 
We state below two localization results which hold in the same regimes. The dynamical local- 
ization property in 7/- is measured in terms of the boundedness in time of all quantum moments 
of the position operator on the lattice. More precisely, for p > we let be the maximal 
multiplication operator such that 



1^ \e^] 



where |j|e denotes the Euclidean norm on 7/. 



for j e TJ^ 



(23) 



For the unitary Anderson model the parameter t takes the role of a disorder parameter. Small 
values of t correspond to large disorder in the sense that is dominated by its diagonal part, 
since Sd{t) tends to the identity as t — > 0. The following result says that in any dimension, 
dynamical localization holds throughout the spectrum of Uu)-, provided t is small enough: 



Theorem 4.1 JE^ [T^/ Consider 
to > such that for all t < to, (y{Uu) 
C < oo such that for all j, k e7/ 



defined by (22), under assumption B. Then, there exists 



a. 



pp 



(Uu,) almost surely. Moreover, there exist a > 0, 



E 



^^P\{ej\Kek)\ 



(24) 



Consequently, for any p > and for any (p in Z (Z ) of compact support, 

sup\\\X\^U^(p\\ < oo almost surely. 



(25) 



Let us consider now the band edge regime. At this point, it is useful to point out that the 
periodicity along the diagonal of the matrix 5* and the definition of make the operator Uuj 
ergodic with respect to the 2-shift in f2 = . By the general theory of ergodic operators, see 
[I3j, it follows that the spectrum of U^^ is almost surely deterministic, i.e. there is a subset S 
of the unit circle such that cr(f4;) = S for almost every u. The same is true for the absolutely 
continuous, singular continuous and pure point parts of the spectrum. Explicitely, there are T,ac, 
Esc and Spp such that almost surely aaciU^)) = Sac, crsc(f4;) = ^sc and cTpp^Uu,) = Epp. Moreover, 
E can be characterized in terms of the support of u and of the spectrum of Sd{t), [19]: 

E = exp {—i supp u) (T{Sd{t)) = {e*" : a G [—dXo, dXo] — supp z/}. (26) 

These facts also hold for the one dimensional half lattice operator (j6]). 

For simplicity, and without loss of generality, we assume that supp u C [— /3, f3] with (3 E (0, vr) 
and —/3, /3 G supp z/. Furthermore, we will work under 
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Assumption C: 

l3 + dXo< TT. (27) 



By (26), this implies tlie existence of a gap in the almost sure spectrum oi U^, 

{e'^ : ^ G {dXo + /3, 27r - rfAo - /?)} n S = 0, (28) 

and that e*^''^''"'"^) and ^^('^^-'^■^0-/3) g^^.^ j^^n^j edges of S. In any dimension, and for any disorder, 
the result below states that localization takes place at the band edges, at arcs denotes by / in 
figure |3} 




Figure 3: The spectrum of and regions of localization /. 

To focus on specific parts of the spectrum of U^, we introduce spectral projectors Pl^^] on 
intervals [a, b] C T, by ^ = E"^([e", e*'']), where dE'^ is the spectral family of and [e'", e'''] 
is a positively oriented arc on §. 

Theorem 4.2 Consider f/^ defined by ^2^ , under assumptions B and C. Then, there exists 
7 > such that for the interval [a, b] = [dXo + /3 — 7, c/Aq + P] it holds: 

(a, 6) n S 7^ and (a, b) n Econt = 0, (29) 

where I^cont = ^sc^^ac- In other words, almost surely P\^a.wP^ pure point spectrum. Moreover, 
there exist constants C < 00 and a > such that 

E[sup|(e,|f/:P[;,]e,)|]<Ce-"l^--^l (30) 

for all j, k And, consequently, for any p > and for any ip in (^{T/) with compact support, 

sup|||X|P[/^Pj'^;,]^|| < 00 almost surely. (31) 

The same is true for the interval [a, b] = [27t — dXo — /3, 27r — dXo — /3 + 7]. 
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5 Quantum Walks in Random Environments 



We now turn to another type of discrete quantum dynamics which can be shown to display 
locahzation due to the presence of disorder. Quantum walks have become a popular topic of 
study due to the role they play in theoretical quantum computing, see e.g. [28], [21], [33], [30]..... 
to their use in the description of effective dynamics of quantum systems, see [1], [22], [IQ], and 
to the nice mathematical properties they have, [2], [2S], [TTj . 

Let us consider the simplest instance of a quantum walk, i.e. a quantum walk on Z. Such 
walks simply describe the discrete dynamics of a quantum particle with spin. In this context, 
the spin state is often called coin state. The Hilbert space is thus 

■H = C^®l\Z). (32) 

We denote by {| t); I i)} ^ canonical basis of and we denote the (position) canonical basis 
vectors of /^(Z) by {|^)}, n G Z. The time-one dynamics of the system is composed of two steps: 
a unitary evolution of the spin alone by means of a unitary operator C on C^, followed by the 
motion of the walker, conditioned on the spin state. More precisely, if the spin is pointing up 
the walker moves to the right one step, and if the spin is pointing down the walker moves to the 
left. The latter step is determined by the action implemented by the unitary operator 

S = J2{Pt^\k + +Pi^\k- l){k\} (33) 



where we have introduced the orthogonal projections 

Pt = |t)(tl and Pi=\l){l\. (34) 

Altogether, the one step dynamics consists in shuffling the spin variable and then performing the 
spin dependent shift S: 



U = S{C^ I) with C 



a b 
c d 



s.t. C* = C-^. (35) 



The evolution operator at time n then reads f/". 

Hence, if one starts form the state |t) ® \k), the (quantum) probability to reach, in one time 
step, the site \k + 1) equals |ap whereas that to reach \k — 1) equals 1 — |ap. Similarly, starting 
from I 4,) ® \k), the probability to reach the site \k + 1) equals |6p and that to reach \k — 1) 
is 1 — The similarity in interpretation with a classical random walk explains why the spin 
variable and the operator C are called the coin states and coin operator. Despite the similarity 
of this dynamics with that of a classical random walk, there is nothing random in the quantum 
dynamical system at hand. The dynamics is invariant under translations on the lattice Z, which 
implies ballistic transport. 

More quantitatively, let X = I ® x denote the operator defined on (its maximal domain in) 
® /^(Z), where x is the position operator given by x\k) = k\k), for all A; G Z. For any p > 0, 
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n G Z, and any with compact support, we consider the quantum mechanical expectation of 
the operator X at time n by 

(XP)^(n) = ((/?|t/-"XP[/"^) = ||XP/2f/"y,||2. (36) 
The analog definition holds for {\X\P)^{n). By Fourier transform methods, one gets 
Lemma 5.1 For any y9 G "H with compact support, 

n->oo n 

with B = iff' C is off' diagonal. 

When C is off diagonal, complete localization takes place. 

A quantum walk in a non-trivial environment is characterized by coin operators that depends 
on the position of the walker: for every A; G Z we have a unitary on C^, and the one step 
dynamics is given by 

U = J2 {P^Ck (^\k + l){k\ + P^Ck ®\k- l){k\} . (37) 
feez 

We consider a random environment in which the coin operator Ck is a random element of U{2), 
satisfying the following requirements: 

Assumption D: 

(a) {Ck}kez are independent and identically distributed f/(2)-valued random variables. 

(b) The quantum amplitudes of the transitions to the right and to the left are independent 
random variables. 

(c) The quantum transition probabilities between neighbouring sites are deterministic and inde- 
pendent of the site. 

There are of course plenty of meaningful ways to define a (random) environment for a quantum 
walk, see e.g. |311 [231 ESI ESI |32]. Assumption D is motivated by simphcity and by proximity 
with the classical random walk. It turns out this choice actually dictates the form of the random 
coin operators as follows. 



Lemma 5.2 IF?] / Under Assumption D, the operator U^^ defined by (37) is unitarily equivalent 
to the one defined by the choice 



where < t, r < 1 and + = i (33) 



and {ujl}kez U {cu^jfcgz are i.i.d. random variables defined as in up to multiplication by 
global deterministic phase. 
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Let be the one step dynamics of a quantum walk in a random environment defined by 
(37) with Cfc, k E Z given by (|38), where {a;*}fcez,#e{t,4} i-i-*^- random variables defined 

in (5]), distributed according to an absolutely continuous measure u on T. Then a statement 
equivalent to Theorem 2.1 in this context holds. 



Theorem 5.3 ^21 J Assume B holds for the distribution dv. Then, for any t G (0, 1), 

(^(UiS) = (^ppiU^) almost surely. 
Moreover, there exist C < oo, a > such that for any j, G Z and any a,T E {f, 1} 



E 



sup \{a®j\U^T®k)\ 



and, for any p > 0, almost surely. 



sup{XTJn) < oo. 



(39) 



(40) 



The similarity in this result and Theorem 2A_ stems from the similarity of the random unitary 
operators in the two cases considered. More specifically, Lemma [5.2| shows that, up to unitary 
equivalence and multiplication by a global phase, f/^ has the following representation in the 
ordered basis {ckjkez = {• • • , 1 1) ® 1^ - !)> I i) ® - 1), 1 1) ® !""■), 11) ® . . .}, 



Uuj = D^S, with S 



t 


r 
t -r 

t 



t 





-r 

t 



t 





(41) 



Here the diagonal of S consists of zeroes and the labeling of the basis is such that the odd rows 
contain r, t and the even rows contain t, — r. Moreover, upon relabeling the indices of the random 
phases, D^^ is diagonal with i.i.d. entries, = diag(. . . , e" 



-10' 



Note that since the random operator at hand differs from that of Remark 2.2 by the form of 



the deterministic matrix S, the localization result stated in Theorem 5.3 requires the separate 
analysis provided in [2T] . 



6 Methods 

Now that we have described several similar random unitary operators appearing in the study of 
different quantum models, we want to address the methods used to derive dynamical localization 



13 



results for these operators. The paper [18] is devoted to a detailed and hopefully pedagogical 
exposition of these methods, so we only point out here the main steps of the analysis. As 
mentioned already, the analysis draws on the similarity of these random unitary models with the 
self-adjoint discrete Anderson model. Actually, our approach to localization proofs will be via 
a unitary version of the fractional moment method, which was initiated as a tool in the theory 
of selfadjoint Anderson models by Aizenman and Molchanov in [7j. Dynamical localization will 
follow as a general consequence of exponential decay of spatial correlations in the fractional 
moments of Green's function. 

Let us consider a random unitary matrix with a band structure in a distinguished basis 
{efejfcezfe of /^(Z'') of the form 

U^ = D^Sd, withZ}^= diag{e-'^^} (42) 

where the random phases {e~*^fe j/.g^'* satisfy assumption B (adapted to the d-dimensional setting) 



and the matrix Sd is a d-fold tensor product of the five-diagonal unitary operators (13) invariant 
under the 2-shift. Again, some results hold under weaker hypotheses, but we stick to our setting 
in order to keep things simple. 
Let 

GUkJ;z) = {ek\{U^-z)-'ei) (43) 

be the Green function of Uuj defined for z G p(US), the resolvent set of ?7^. Now, the structure 
of is such that a modification in one of the random parameters corresponds to a rank one 
perturbation of the original operator. This leads to the observation that while the Green func- 
tion becomes singular as z approaches the spectrum of f/j^, these singularities are fractionally 
integrable with respect to the random parameters: for s G (0, 1) the fractional moments of the 
resolvent^ E(|G'(/i;, /; z)Y), have bounds which are uniform for z arbitrarily close to the spectrum. 
This is the content of our first result. 

Theorem 6.1 Suppose assumption B holds for the random variables {0k}kei'^- Then for every 
s G (0, 1) there exists C{s) < oo such that 




\G^{k,l;z)\'diy{ek)diy{ei)<C{s) (44) 

for all 2; G C, \z\ 7^ 1, all k,l G Z'^, and arbitrary values of 6j, j ^ {k,l}. Consequently, 

E{\G^{k,l;z)n<C{s), (45) 

for all z ^ C, \z\ ^ 1. 

Remark 6.2 The proof of this general result makes use of the fact that the measure du has a 
density in . 

Then, the goal is to make use of the specificities of the model under study to identify regimes 
or situations where the fractional moments E(|G(/c, /; 2;)|*) are not just uniformly bounded, but 
decay exponentially in the distance between k and /. The following general result shows that this 
can be used as a criterion for dynamical localization of U^. 
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Theorem 6.3 Suppose assumption B holds for the random variables {9k}k&'^ ^^^'^ that for some 
s G (0, 1), C < oo, a > 0, £ > and an interval [a, h] G T, 

E(|G'(A;,/;;z)n < Ce-"l^-^l (46) 

for all k,l El/ and all z E C such that 1 — 5 < |2;| < 1 and argz G [a,b]. 
Then there exists C such that 

E[ sup |(e,|/(f/jPf^,,]Q)|] <C'e-"l'=-'l/^ (47) 

/eC(S) 

11/11 oo<l 

for all kj eZ'^. 



Remark 6.4 That the estimate (^7) implies almost sure spectral localization on {a,b) can be 
shown by means of arguments of Enss-Veselic ITB^ on the geometric characterization of bound 
states. Also, (4^) directly prevents the spreading of the wave function over all times, in the sense 
that for all p > 0, sup^g^ || iXl^t/^Pj*^ < oo almost surely. Both these facts are explicitly 
shown in 

Note that speciahzing to the case f{z) = z"-, with n G Z, we get the locahzation resuhs stated 
in the previous sections. 



The proof of Theorem 6.3 requires a hnk between the fractional powers of the resolvent and 
the resolvent itself, so that some functional calculus can be applied to control operators of the 
form f{U), for certain continuous functions / : § — t- C. This is done in two steps. The first one 
is an estimate on the expectation of the square of the Green function in terms of the expectation 
of fractional powers of the Green function. This step is equivalent in our unitary framework to 
the second moment estimate proven by Graf in (TB] for the self-adjoint case. 

Proposition 6.5 Assume B. Then for every s G (0, 1) there exists C{s) < oo such that 

E{{1 -\znG{k,l;z)\')<C{s) J2 n\G{m,l;z)n (48) 

\m-k\<A 

for all \z\ < 1 and k,l E Z'^. 

Remark 6.6 The fact that the sum in the right hand side of the inequality only carries over 
indices m a finite distance away from k is a direct consequence of fact that the deterministic 
operator S has a band structure. 

The second step consists in reducing bounds for f{U) to bounds on resolvents by means of 
the following result. 



Lemma 6.7 



1 - /■^'^ 

f(U) =w - lim / (U - re'^)-\U-^ - re^'^Y^ f{e'^)de (49) 



for f E C{S) and U a unitary operator. 
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Remark 6.8 This formula is a consequence of the representation of non-negative Borel measures 
on T hy Poisson integrals. This can he seen by considering the non negative spectral measure 
dfi^p on the torus T associated with a normalized ip ^Ti such that {(f\U(p) = J^e^^dfi^pla), and 

(1 - r'){v\{U - re^')-\U-' - re-^)-V) = / , (50) 

JT F — re I 

For any f E C{S), we thus have 

r^TT r 1 2 Jf) 

WiUM^m^ l X^.j37^*.W/(e-% (51) 

and one concludes by polarization. 



If the fractional moments of the resolvent are exponentially decaying, i.e. if (46) holds, so 



is the left hand side of (48). Then, considering matrix elements of (49) and applying Fatou's 



lemma and Cauchy Schwarz, one derives the upper bound (47), as shown in 



We have seen that showing dynamical localization for a concrete model amounts to proving 
that the fractional moments of the resolvent are exponentially decaying, i.e. that (46) holds. 



This has been done in different ways for the different regimes and models considered. We shall 
not attempt to explain in details how of this technical task is achieved in the models considered 
above, but we just want to describe the methods employed to do so. 

For one dimensional models, either on /^(N) or on /^(Z), one studies the generalized eigen- 
vectors of the problem, i.e. the solutions to U^ji' = zip in l(N) or /(Z). Because of the band 
structure of the operator [/^, the generalized eigenvectors are obtained by means of a transfer 
matrix formalism and their behavior at infinity is controlled by the associated Lyapunov ex- 
ponent. Exploiting the way the randomness appears in the model, one then shows that the 
Lyapunov exponent is positive and continuous in the spectral parameter z, in a neighborhood 
of the unit circle. Then, by making use of the expression of the Green function in terms of 



certain generalized eigenvectors, one shows that (46) holds throughout the spectrum, and for all 
values of the parameter t G (0, 1). This strategy was implemented in [18] for the magnetic ring 
model and for the one dimensional unitary Anderson model, and in for the quantum walks 
in random environments models. Previous studies of the properties of the Lyapunov exponents 
for these models were performed [HllinillZI, which lead to spectral localization results by spectral 
averaging, according to a unitary version of the argument of Simon- Wolff, |36j . 

For the (i-dimensional unitary Anderson model, the large disorder regime was addressed in 
[20] . It was shown in this paper that estimate (46) holds in any dimension, provided t is small 
enough. To prove this estimate, the similarity in the way the randomness appears in the model 



(22) with the discrete Anderson model was used explicitly. The analysis is based on estimates on 



the expectation of the resolvent equation raised to a fractional power s, on rank one perturbation 
formulas and on a so called "decoupling Lemma", similar to the one shown in [I] for the self- 
adjoint case. This leads to an inequality satisfied by the function < f{k) = K{\F{k,j; z)\'^) in 
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where F{z) = Uuj{Uuj — z)~^ = I + z{Ui^ — z)~^ is essentially equivalent to the resolvent. 
This inequality says that f{k) is smaller than a 2;-independent constant times the weighted 
average of its values around k, with weights given by the matrix elements of Sd- The structure 
of Sd and dependence in t of its matrix elements then imply the sought for bound, for t small 
enough. 

The band edge regime for the dimensional unitary Anderson model was tackled in the paper 
[T8] . adapting the general strategy provided in [6]. This regime, which is the most challenging 
to cover, requires getting finite volume estimates on the resolvent, close to the band edges. A 
first step consists in defining the restriction [/a;|A(L) of to finite boxes A(L) C Z'^ of side 
length L by means of appropriate boundary conditions which make this restriction unitary and 
imply certain monotony properties of the spectrum as boxes are spit by adding more boundary 
conditions. Then, one needs to get accurate probabilistic bounds on the size of the resolvent of 
this restriction, when the spectral parameter z is close to the band edges. It requires showing 
that when L becomes large, the probability to have eigenvalues a distance smaller than 
away from the band edges is of order e~'^^°' , for < /3 < 1 and a,7 > 0, i.e. a Lifshitz tail type 
estimate. Then a decoupling lemma with an iterative argument allows us to prove the bound 



(46) for the infinite volume operator f/^, in a non-empty neighborhood of the band edges. 

Finally, we would like to mention that there is at least one more popular model in condensed 
matter physics whose dynamics reduces to the study of a discrete time quantum dynamics gen- 
erated by a random unitary operator with a band structure: the Chalker Coddington model 
and its variants, see [lO]. This model can be thought of as a unitary equivalent of the discrete 
Schrodinger equation on a finite width two-dimensional strip. Some progress was made recently 
about the properties of this model in [1] . But the focus of this work is more on the analysis of 
the associated set of Lyapunov exponents than on dynamical localization aspects. This is why 
we didn't provide a description of the Chalker Coddington model in these notes, eventhough it 
certainly belongs to the family of unitary random operators presented here. 

Acknowledgements: It is a pleasure to thank Bob Sims and Daniel Ueltschi for the 
invitation to the perfectly organized "Arizona School of Analysis with Applications 2010", where 
part of this material was presented. 
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